Electronic Structure of Strongly Correlated Systems Emerging from Combining 
Path-Integral Renormalization Group with Density Functional Approach 
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A new scheme of first-principles computation for strongly correlated electron systems is proposed. 
This scheme starts from the local-density approximation (LDA) at high-energy band structure, while 
the low-energy effective Hamiltonian is constructed by a downfolding procedure using combinations 
of the constrained LDA and the GW method. Thus obtained low-energy Hamiltonian is solved 
by the path-integral renormalization-group method, where spatial and dynamical fluctuations are 
fully considered. An application to Sr2VC>4 shows that the scheme is powerful in agreement with 
experimental results. It further predicts a nontrivial orbital-stripe order. 
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First principles methods for electronic structure cal- 
culations have been extensively applied to various sys- 
tems in the last few decades. The density functional 
theory (DFT) supplemented with the local-density- 
approximation (LDA) 0, 0] is one of the most successful 
schemes. However, in strongly correlated electron sys- 
tems such as transition metal oxides and organic com- 
pounds, LDA fails even at qualitative levels of simple 
systems and remain challenges in terms of first-principles 
calculations. The failure is ascribed to the correlation ef- 
fects generating large spatial and dynamical fluctuations 
beyond the applicability of existing schemes. 

In electronic structure calculations, DFT and wave- 
function methods are two typical approaches 3]. The 
rapid development and recent extensive applications of 
DFT approach rely on its less computation time while 
a systematic improvement of its accuracy is not an easy 
task. On the other hand, the wavefunction methods in- 
cluding the Hartree-Fock and variational schemes cost 
much longer computational time, while it offers a better 
accuracy particularly for strongly correlated systems as 
we will discuss later. 

After considering this contrast, an optimal choice of 
the algorithm would be a hybrid method. Since the elec- 
trons far from the Fermi level do not show serious fluctu- 
ation effects, LDA and GW HG| offer a reasonable and 
efficient computational scheme there. If the tracing out 
and the elimination of the high-energy degrees of free- 
dom by a downfolding scheme would result in the effec- 
tive low-energy Lagrangian or Hamiltonian, they can be 
treated by an accurate wavefunction method. Such suc- 
cessive eliminations of the high-energy part have a con- 
ceptual similarity to the renormalization group method. 
In this paper, we propose a downfolding scheme to derive 
the low-energy effective models from the higher-energy 
structure, together with a reliable low-energy solver for 
the downfolded effective model, as a single framework 
of a first-principles method. This offers a powerful first- 
principle method for strongly correlated electron systems. 
To show the performance, we apply the present scheme to 
Sr2VC>4, which reveals intriguing properties of this com- 
pound. 

Since the electrons have kinetic and interaction ener- 



gies, the downfolding procedure basically consists of two 
parts. One is to derive the screened Coulomb interaction 
for low-energy electrons after eliminating the high-energy 
degrees of freedom. The other is to take into account the 
effect on the single-particle kinetic part, which can be 
expressed as the self-energy effect from the high-energy 
electrons, if the quasi-particle description is justified. 

After eliminating the high-energy degrees of freedom, 
the effective model in general has a frequency dependence 
represented by a Lagrangian. However, in the low-energy 
region, the Hamiltonian approach by replacing the dy- 
namical Coulomb interaction with the static one still of- 
fers an efficient and essentially correct framework if a 
certain condition is satisfied 

Recently, the dynamical mean-field theory (DMFT) || 
was employed as a low-energy solver for finite tempera- 
ture properties 0- Since the DMFT ignores the spatial 
correlations, its refinements by using the cluster or cel- 
lular algorithms [Tol ITU were attempted. In this paper, 
we alternatively solve the low-energy effective Hamilto- 
nian by the path-integral renormalization group (PIRG) 
method |l2j implemented with quantum number projec- 
tion algorithm 13] . The PIRG scheme allows treating 
equally the spatial and dynamical fluctuations in a con- 
trollable way and hence makes it possible to obtain the 
ground state of the effective Hamiltonian exactly in a 
numerical optimization scheme. 

We now start with the downfolding procedure. We first 
perform the LDA calculations based on the LMTO basis 
functions ^4(. In many correlated electron systems, the 
bands close to the Fermi level are relatively well sepa- 
rated from the high-energy band structure as in transi- 
tion metal oxides, where the bands near the Fermi level 
consist of 3d character mainly derived from transition- 
metal atomic orbitals (and frequently hybridized with 
oxygen 2p atomic orbitals). In this circumstance, the 
whole bands may be downfolded to an effective model 
consisting only of the transition metal 3d bands (or ad- 
ditionally with the oxygen 2p bands). If the crystal field 
splitting is strong, it may be downfolded even to a set of 
the bands closest to the Fermi level, for example, the ti g 
or e g bands only, in a cubic environment. Here we report 
the first attempt of the downfolding along this line. 
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FIG. 1: (color) (Right panel) :Comparison of LDA dispersions 
for Sr2VC>4 computed from the LMTO basis functions (dotted 
curve) with the downfolded tight-binding bands (dot-dashed 
curve). (Left panel) :The corresponding LDA density of states. 

When we restrict the degrees of freedom to the isolated 
LDA bands closest to the Fermi level, we first need to 
reduce the Hilbert space. With the notation for the states 
having the maximal weight near the Fermi level as {\d)} 
and the rest of the basis functions as {\r)}, the whole 
Hilbert space may be spanned as {|x}}— {|^)}©{l r )}- A 
typical example of Si'2 VO4 is shown in Fig.^ In this case, 
three isolated bands closest to the Fermi level mainly 
consist of the 3d ti g states, which defines the Hilbert 
subspace {\d)} in our downfolding procedure. The LDA 
eigenvalue equation can be rewritten as 

(H dd - u)\d) + H dr \r) = 0, (1) 
H rd \d) + (H rr -u)\r) = 0. (2) 

We eliminate the subspace \r), which yields the effective 
w-dependent Hamiltonian for only d-subspace: H d ^{u) = 
fjdd _ H dr (H rr - uj)- 1 ^. From the overlap matrix 
S{uj) = 1 + H dr (H rr - cj)- 2 H rd , with (d\S\d)=l, the 
tight-binding effective Hamiltonian, h is obtained after 
the orthonormalization of the vectors \d)— ► |d)=S 1 / 2 |d) 
and fixing the energy u in the center of gravity of the 
band of our interest: h = S~ 1/2 H d ^S^ 1/2 . After Fourier 
transformation to the real space, h defines the Wannier 
basis. The band structure obtained after the elimination 
of the higher-energy bands shows an excellent agreement 
with the LDA band computed from the LMTO basis as 
we see in Fig. ^ This formalism has generality and can 
be easily extended to a more complex band structure |15| . 

We now derive an effective Coulomb interaction W r 
among electrons at the isolated t2 g Wannier orbitals af- 
ter considering the screening by other bands. In princi- 
ple, W r may be derived by a full GW scheme 0. How- 
ever, for practical use, the whole GW calculation re- 
quires large computation time. Since the screening from 
the polarizations of high-energy electrons does not have 
strong fluctuations, it may alternatively be replaced by 
the constrained-LDA (C-LDA) method following the pro- 
cedure in the literature |l6|. Then we take two steps 
in deriving the screening. First, we compute the in- 
teraction W r i which takes into account the screening 



by the electrons residing at atomic orbitals other than 
3d. This step is performed by the C-LDA method [l6| . 
When the higher-energy bands are well separated, we ex- 
pect that the frequency dependence is small around the 
Fermi level and it may well be calculated by the C-LDA 
scheme by starting with the basic definition of the 
effective Coulomb interaction formulated by Herring |17| : 
Unn> = E[n R + l,n R > -1] -E[n R , nn>], which is nothing 
but the energy cost for moving a 3d-electron between two 
atoms, located at R and R', and initially populated by 
7iR=nR,'=n electrons. Note that at this stage, the screen- 
ings of 3d-3d interactions by the electrons at the same 
3d atomic orbitals are excluded, because the 3d atomic 
orbitals do not hybridize each other in C-LDA scheme. 
For Sr2V04, the onsite Coulomb interations derived from 
[/rr/ is C/~11.3 eV. For comparison, the bare Coulomb 
interaction is 21.8 eV. Thus, the screening of 3d-3d in- 
teractions by non-3d electrons is strong. In addition to 
the parameter U, the intra-atomic exchange coupling J 
can be easily calculated in the framework of C-LDA ap- 
proach. By using these two parameters, one can restore 
the full matrix W T \ of screened Coulomb interactions p^|. 

It is known that LDA does not take into account the 
self-energy effect and fails in reproducing a gap structure 
at the Fermi level [5j- LDA also ignores energy shifts in 
core electron levels arising from the self-interaction ef- 
fect. However, in the present scheme, the bands near 
the Fermi level are left for a more refined PIRG method, 
while the core electrons hardly polarize and does not con- 
tribute to the screening. The valence electrons give the 
major contribution to the screening in this C-LDA treat- 
ment and their energies are well described by LDA with 
few self-energy corrections in many cases including tran- 
sition metal oxides Therefore, although C-LDA is 
not complete, in the practical usage, it offers a reasonable 
and efficient way of counting the screening effects by the 
valence electrons. 

In the second step, we use the GW scheme @ by tak- 
ing W r i as if it is the starting bare Coulomb interaction. 
Then the RPA screening produces W r (uj) = W r \/{\ — 
PdrW r {). The whole polarization of the 3d atomic or- 
bitals in the LMTO basis is given by Pd = Pdr + Pt2g, 
where Ptig is the polarization purely from the tig LDA 
band, while P ( i r represents the rest of 3d atomic or- 
bitals contribution contained in the e g LDA band as well 
as in the component hybridized in the oxygen 2p LDA 
band. Since the identity W[u) = W rl /(1 - P d W r i) = 
W r /(1 — Pt2gW r ) holds, W r plays the role of the effec- 
tive interaction in the downfolded Hilbert space of the 
t2g LDA band. We also note that some of the oxygen 2p 
LMTO component is contained in the 3d ti g Wannier or- 
bital. The screening by these components is already taken 
into account by the C-LDA scheme, while it should be 
excluded in the low-energy model for the 3d tig Wannier 
orbitals. In practice, however, this contribution to the 
screening turns out to be negligible. Although the whole 
GW scheme considering P t 2 g results in W, we consider 
the interaction within the <2g LDA band more accurately 
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FIG. 2: (color) Diagonal matrix elements of the screened 
Coulomb interaction W r (ui) (left panel) and the self-energy, 
(right panel). The inset shows the high-frequency part 
of W r (uj) (the lines corresponding to the xy and yz orbitals 
are indistiguishably close). The real and imaginary parts are 
shown by solid and dashed curves, respectively. 



by the low-energy solver as we explained above. Then in 
the effective low-energy model for the basis represention 
of the t2 g Wannier orbital, the screened Coulomb inter- 
action is given by W r . In general, W r becomes frequency 
dependent (Fig. 0). If the frequency dependence is small 
within the range of the ti g bandwidth, we are allowed to 
take the low-frequency limit W r (u) = 0) as the interac- 
tion part U in the low-energy effective Hamiltonian [(|. 

In the downfolding process, the kinetic-energy part is 
modified through the self-energy £(fc,o>), which can be 
evaluated in the GW approximation Here the self- 
energy effect from 3c? and 2p atomic orbitals is consid- 
ered while the polarization of the higher energy bands 
is small and we neglect it. We also take into account 
the self-energy arising from the dynamical part of the 
Coulomb interaction W r (ui) — W r (ui = 0) through the 
GW scheme l(i . Such a self-energy effect coming from 
the higher energy part of W r mainly appears through 
ReE, and contributes to the renormalization factor Z = 
(l-5E/9w)-I (~ 0.8 for Sr 2 V0 4 ). The low-energy part 
of ImE is small and can be ignored (Fig[3J). Then the ef- 
fective Hamiltonian is reduced to a multi-band Hubbard 
model in the Wannier representation for the ti g band: 



H = 



m,m f ,£T 



'13 



i,ct,{3 
7,5 



(3) 

where c| 7no .(c, m(T ) creates (annihilates) an electron with 
the spin a = (f, J,) at the t2 g Wannier orbital m = 
(xy, yz, zx) of the site i and n imcr = c\ ma c tm(J . The Greek 
symbols stand for the combination (rn, a) of the indices. 

We now discuss the PIRG method as a solver of the 
low-energy effective Hamiltonian. The basic method |l2^ 
is to perform the path-integral operation by following 
the principle that lim p ^ 00 [exp(— tH)} p \&) generates the 
ground state, with a proper starting trial wavefunction 
$). Through the summation over the Stratonovich vari- 
ables in the Stratonovich-Hubbard transformation of the 
interaction term, the obtained state is given by a linear 



combination of increased number of the nonorthogonal 
Slater-determinant basis functions. The optimized ba- 
sis construction with systematic extrapolation to the full 
Hilbert space leads to an accurate estimate of the ground 
state. It is also crucial to implement the quantum num- 
ber projection in restoring the symmetry of the Hamil- 
tonian and also in improving the accuracy fl3T ]. A typi- 
cal accuracy is seen in the ground-state energy estimate 
of —1.85790 ± 0.00002 in comparison with the quantum 
Monte Carlo result of -1.8574 ± 0.0014 for the standard 
Hubbard model at half filling on a 6 x 6 square lattice 
with the nearest neighbor transfer t = 1 at the onsite 
interaction [7 = 4 \l3\ . 

The present first-principles method combined with the 
PIRG solver is now applied to calculate the electronic 
structure of Sr2V04. This compound has a layered per- 
ovskite structure and bears a intriguing character be- 
cause this has one 3d electron per V site (d 1 system) 
with strong two-dimensional anisotropy and has a dual 
relation to the one 3d hole per Cu sites (d 9 system) in the 
mother compounds of the cuprate superconductors. The 
duality is not complete because, in Sr2V04, the orbital 
degeneracy of d 1 electron perfectly remains between d yz 
and d zx orbitals. The crystal field splitting of d xy orbital 
is also rather small (~ 0.08 eV in the LDA calculation). 

In case of Sr 2 V04, after the downfolding, the onsite 
interactions among the Wannier orbitals of intraorbital 
xy, yz(zx) and interorbital xy-yz(xy-zx) and yz-zx com- 
binations are U = 2.77,2.58,1.35 and 1.28 eV, respec- 
tively. The onsite exchange interactions between xy- 
yz(xy-zx) and yz-zx orbitals are 0.65 and 0.64, respec- 
tively. The transfer between the nearest neighbor xy- 
xy,yz-yz and zx-zx orbitals in x direction are -0.22, - 
0.05 and -0.19 eV, respectively. Although the transfer 
between the next-nearest-neighbor pairs are considerably 
smaller, they have also been taken into account in our 
calculations. In order to monitor the effect of Coulomb 
interaction, we examine the scale-factor dependence by 
multiplying all the matrix elements U a f3 7 s in Eq.(P]) with 
a factor A. Namely, the realistic value in the downfoldcd 
Hamiltonian corresponds to A = 1. 

Sr2V04 was first experimentally studied by Cyrot et al. 
and Nozaki et al. |20ll2l| and recently thin film formed on 
the substrate LaA103 was studied by Matsuno et al. [2^] • 
Transport and optical properties indicate that this com- 
pound has either a very small Mott insulating gap or 
semiconducting property with rapidly increasing resistiv- 
ity with decreasing temperature. The gap amplitude is 
nearly zero and it can be easily metallized by Sr dop- 
ing |22| . The magnetic susceptibility appears to show 
presumable antiferromagnetic transition at around 50K. 

In sharp contrast to this nearly insulating transport 
properties, the LDA calculation predicts a good metal- 
lic behavior (Fig^J. On the other hand, the Hartree 
Fock calculation predicts a clear ferromagnetic insulat- 
ing phase at the realistic parameter values (see Fig. |3J). 

The PIRG results indicate that after quantum fluctua- 
tions taken into account, the ground state becomes close 
to the metal-insulator transition point and becomes to- 
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FIG. 3: Lowest energies per unit cell of total S=0 (open sym- 
bols) and ferromagnetic states (filled symbols) by quantum- 
number projected Hartree-Fock (triangles) and PIRG (circles) 
calculations for the downfolded model of S^VCU. The en- 
ergy unit is eV. (Inset) :Lowest energies per unit cell of metal- 
lic (open squares) and insulating states (filled squares) by 
quantum-number projected PIRG calculations for the down- 
folded model of Sr 2 V0 4 . 

tal spin S = in contrast to the Hartree-Fock predic- 
tion as we see in Fig. [21 The Mott transition occurs at 
A ~ 0.95 as indicated by the jump of the double occu- 
pancy D and the realistic parameter is on the verge of 
the first-order Mott transition. The double occupancy is 
defined by D = Q2 m<7 n.i ma -). This is consistent with the 
experimental result. Furthermore, the spin and orbital 
are strongly coupled each other and the PIRG results in 
the spin-orbital order shown in Fig. 21 with a remarkable 
and nontrivial orbital-stripe structure for the insulating 
phase, while it remains paramagnetic in the metallic side. 
The realistic value of the relativistic spin-orbit coupling 
does not seem to alter the present conclusion. Several 
possible interlayer configurations are degenerate within 
the accuracy (~ 0.005eV) of our calculation. At least, 
the antiferromagnetic state is consistent with the experi- 
ments with insulating behavior, while the configurations 



of spin-orbital order are not experimentally available so 
far and is desired to test whether this prediction is correct 
or not. More detailed theoretical study will be reported 
elsewhere p3| . 

It has turned out that this compound shows very se- 
vere competitions. First, it lies on the verge of the Mott 
transition. Second, the ferromagnetic state is rather close 
in energy to the true ground state with the antiferromag- 
netic order. Third, candidates of configurations for the 
spin-orbital order are in severe competitions each other 
in the order of 100K in energy. The available experi- 
mental results are consistent with our present results in 
contrast to the LDA and Hartree-Fock results. Clearly 
we need more analyses on this compound. However, 
all of the above indicate that our present approach of 
PIRG combined with the downfolding by using the LDA- 
GW scheme offers a promising computational method for 




FIG. 4: Ordered spin-orbital patterns in plane for Sr2V04 in 
the Mott insulating ground state with the unit cell of 2 x 4. 
Ordered spin moment is proportional to the length of arrows. 



strongly correlated electron systems. From the viewpoint 
of the computational methods, this compound appears to 
offer a very severe and good benchmark test for the ac- 
curacy in taking account of the correlation effects. 

The authors are grateful to F. Aryasetiawan for fruitful 
discussions and also thank Y. Tokura for clarifications of 
experimental results on Sr2V04. 
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